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Abundant H2D + in the pre-stellar core L1544 
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Abstract. We have detected the 372 GHz line of ort/io-rFD" 1 " towards the pre-stellar core L1544. Fhe strongest 
emission (T m b ~ 1 K) occurs at the peak of the millimeter continuum emission, while measurements at offset 
positions indicate that Bi2D + is confined within ~ 20", where CO is highly depleted. The derived H2D + abundance 
of ~10~ 9 is comparable with previous estimates of the electron abundance in the core, which suggests that H2D + 
is the main molecular ion in the central 20" (2800 AU) of L1544. This confirms the expectations that H2D + is 
dramatically enhanced in gas depleted of molecules other than H2. The measured abundance even exceeds the 
present model predictions by about a factor ten. One possibility is that all CNO-bearing neutral species, including 
atomic oxygen, are almost completely (<; 98%) frozen within a radius of ~2800 AU. 
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1. Introduction 

The last few years have seen a boom of the studies of 
molecular deuteration in star forming regions, triggered 
by the discovery of a large fraction of doubly deuter- 
ated formaldehyde in the low mass protostar IRAS16293- 
2422 HCeccarelli et al. 19981 ILoinard et al. 2000|l . where 
the measured D2CO/H2CO ratio is 25 times larger than 
in Orion IjTurner 1990jl . It is now clear that IRAS16293- 
2422 is not unusual, for all Class protostars show sim- 
ilar or larger D 2 CO/H 2 CO ratios IjLoinard et al. 20 02). 
Formaldehyde is not unusual either, as other molecules 
present astonishing deuteration degrees. Notable ex- 
amples are methanol, whose deuterated forms are as 
abundant as the main isotopomer (Pari se~et al. 20 02). 
or ammonia, where not only the doubly deuter- 
ated form URoueff et al. 20001 ILoinard et al. 2001|l, but 
also the triply deuterated form Ijvan der Tak et af~2 002 
ILis et al. 2002J) has been detected. 

The picture emerging from those studies is that 
this extreme molecular deuteration starts during the 
pre-collapse phase, when the gas is cold and dense 
(Ceccarel li et al. 2 001). Whether formed in the gas 
phase (e.g. IRoberts fc Millar 2000|l or on the grain sur- 
faces fTiel ens 19831 ICaselli et al. 2002a(l , deuterated iso- 
topomers of neutral species are stored into the grain 



mantles, and released back into the gas phase during 
the collapse phase, as the dust is heated by the new 
born star. This picture has been substantially confirmed 
by the recent observations of Bacmann et al. (2003)[ 
who measured D 2 CO/H 2 CO~ 10% in a sample of pre- 
stellar cores presenting large degrees of CO depletion 
HBacmann et al. 2002|l . A decisive factor for the large ob- 
served deuteration seems in fact to be the CO depletion: 
the larger the CO depletion, the larger the molecular 
deuteration. This suggests that the enhanced molecular 
deuteration is the product of gas phase chemistry, given 
that CO depletion leads to an increase of the H 2 D + /H^ 
abundance ratio HDalgarno fc Lepp 1984| ) and thus stim- 
ulates the formation of deuterated molecules, mainly pro- 
duced in ion-neutral reactions with H 2 D + . 

Evidently, the key molecule for fully understand- 
ing the molecular deuteration process is indeed 
H 2 D + . Unfortunately, the only H 2 D + transition 
observable from ground based telescopes, the ortho- 
H 2 D + lio — In transition, lies at 372 GHz, next 
to an atmospheric band which makes the observa- 
tion rather difficult. As a result, the search of H 2 D + 
has been a very frustrating business for about two 
decades (e.g. |Phillips et al. 1985| |Pagani et al. 1992| 
Ivan Dishoeck et al. 19921 IBoreiko fe Betz 1993j) . and 
although several massive protostars have been targetted, 
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H 2 D + has been so far detected, with a relatively low sig- 
nal (main beam temperature, T m b < 0.1 K), in only two 
low mass protostars: NGC1333-IRAS4 IjStark et al. 19990 
and IRAS16293-2422 (Stark et al., in preparation). We 
decided to re-start the H 2 D + search in sources with 
the largest expected H 2 D + /H;j~ ratio, even where the 
excitation conditions of the 372 GHz line may not be 
the most favorable. And, since molecular deuteration is 
indeed dramatically enhanced during the pre-collapse 
phase of low mass protostars, we decided to search for 
H 2 D+ in the best studied pre-stellar core. 

L1544 consists of a dense core surrounded by a 
low density envelope which is undergoing extended in- 
fall IjTafalla et al. 19981 IWTlliams et al. 19990 . From dust 
emission and absorption observations, we know that the 
core has a central density of about 10 6 cm -3 inside a ra- 
dius rfl a t = 2500 AU (the "flattened" region), followed 
by 1/r 2 density fall-off until a radius of about 10000 
AU dWard-Thompson et al.~1999| IBacmann et al. 20000 . 
Molecular species such as CO, CS, and CCS are highly 



Table 1. Results of Gaussian fits to the five H 2 D + spectra 



depleted at densities above 



10 5 



or inside 



radii of ~6000 AU IjCaselli et al. 19991 IQnashi et al. 1 999 
Taf alla et al.l 2002). On the other hand, molecules such 
as NH 3 , N 2 H + and N 2 D + do not show signs of deple- 
tion, probably because of the volatility of their precur- 
sor N 2 flBergin fc Langer 1997| ICaselli et al. 2002bl here- 
after CWZ). Indeed, N 2 H+ and N 2 D+ have been used to 
study gas kinematics (|Caselli et al. 2OO2c0 and the ion- 
ization degree (CWZ) of the high density core nucleus. 
The iY(N 2 D + )/7V(N 2 H+) column density ratio toward the 
dust peak was determined to be ~0.2, in agreement with 
chemical models which take into account differential deple- 
tion of molecular species in centrally concentrated clouds 
IjAikawa et al.l 2001; CWZ). These models predict frac- 
tional abundances of H 2 D + to be ^ 10~ 10 , a value which 
can be tested by the present observations. 

2. Observations and Results 

During two nights in October 2002, we observed the 
dense core L1544 in the lio — In transition of H 2 D + 
(u = 372.42134GHz), using the Caltech Submillimeter 
Observatory on Mauna Kea. The spectra were taken in 
wobbler switching mode, with a chop throw of 300". The 
central position (the peak of the 1.3 mm continuum dust 
emission map of [Ward-Thompson et al. 19990 , as well as 
the 4 positions offset by ±20" in Right Ascension and 
Declination were observed. The rms reached were of the 
order of 110 mK for the central position and 60-90 mK for 
the others. We used as backend an acousto-optical spec- 
trometer with 50 MHz bandwidth. The velocity resolution, 
as measured from a frequency comb scan, is 0.101 kms -1 . 
The beam efficiency at v = 372 GHz was measured on 
Saturn during the observations, and found to be «0.65 
for a 22" FWHM beam size. Pointing was monitored ev- 
ery couple of hours and found to be better than 3". 

In Table \I]we report the integrated intensity, Vlsr ve- 
locity, and line width of the five observed H 2 D+ spectra. 
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In Fig. ^ the peak spectrum is shown together with the 
average spectrum obtained by summing the four off-peak 
spectra. First, the most striking feature of the H 2 D + (lio- 
ln) line at the dust peak is its intensity (~1 K), which 
is at least ten times stronger than any other previously 
detected H 2 D + line ( Stark et al. 1999), suggesting that 
pre-stellar core conditions favour the production of H 2 D + . 
The second important result is the detection of a strong 
signal (~0.4 K) about 20" from the center. Assuming 
that the H 2 D + emission is constant within a certain ra- 
dius r(H 2 D+) and then drops to zero, one can estimate 
r(H 2 D + ) by convolving the predicted H 2 D + integrated 
intensity profile with the CSO beam and comparing it 
with the data. This comparison, shown in Fig. (thin and 
dashed curves), gives r(H 2 D + ) ~ 15"- 20 ", consistent 
with rflat = 17", the size of the "flattened" region ob- 
served in the millimeter continuum emission. In the case of 
a power law distribution of the H 2 D + integrated intensity 
(J T m bdv oc r^ 1 ; dotted curve in Fig. the estimated 
(convolved) size is r(H 2 D+) = 17". This result suggests 
that the bulk of the H 2 D + emission is concentrated within 
the highly depleted high density core nucleus. 

The central line profile shows a dip at the rest ve- 
locity which may be due to the kinematics of the emit- 
ting gas (in analogy with other molecular ion profiles; see 
ICaselli et al. 200"2c)l or it may indicate that the line is op- 
tically thick. We will get back on this issue in a future 
paper where the velocity structure of the L1544 core will 
be included in a detailed radiative transfer code (van der 
Tak et al., in preparation). 

3. Column density and abundance 

As shown in Fig. 2, the H 2 D + emission appears to be 
concentrated within a radius of ~2800 AU. In this re- 
gion, the gas temperature is predicted to be ~ 7 K 
l|Evans et al. 20011 IGalli et al.l 2002), altho ugh molecular 
data suggest larger temperatures (~ 9-10 K; Taf alla et al.l 
2002. IBacmann et al. 2 002). However, because of the large 
amount of gas phase depletion (see Sect.0J, it is likely that 
molecular tracers do not provide information on the gas 
temperature in the central 2800 AU. Assuming LTE con- 
ditions (T cx = 7 K), we can immediately determine the 
optical depth (r) of the H 2 D + (lio~lii) line: 
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Fig. 1. (Top) The H 2 D+(li -ln) line at the dust peak 
of L1544 (RA(1950)=05:01:13.1, Dec(1950)=25:06:35.0). 
The black curve is the Gaussian fit (see Tabled . (Bottom) 
The H 2 D + (li -ln) spectrum averaged in the four posi- 
tions 20" off the dust peak. 
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Fig. 2. Integrated intensity of the H2D + (lio-ln) line as a 
function of projected distance b from the core center. The 
solid and dashed curves show two "step" functions (null 
H 2 D + integrated intensity at impact parameters larger 
than 15 (solid) and 20 (dashed) arcsec) convolved with 
the CSO beam. The dotted curve is a (convolved) power- 



law distribution in the H 2 D + integrated intensity (oc r 
and its size at half-maximum (FWHM) is 34". 
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where T mb =0.91±0.11 K at the peak, T cb = 2.7 
is the cosmic background temperature, and J„(T) = 
T /(exp(T /T) - 1) with T = hv/k = 17.9 K. Assuming 
a spontaneous transition coefficient A — 1.08xl0 -4 s _1 , 
the total column density of ortho-H2D + is: 



so that at the dust peak, where r ~1 and the linewidth is 
0.47 kms _1 (Table[T|), we obtain 2V ort ho = 2.4xl0 13 cm" 2 . 
At the low assumed temperatures, where the ortho to para 
ratio is close to 1 l|Gerlich et al.1 2002). we deduce a total 
H 2 D + column density of ~4.8xl0 13 cm" 2 , about three 
times larger than 7V(N 2 H+) ( see CWZ). 

To gauge the H 2 D + abundance, the above estimate 
of the total H 2 D + column density has been divided by 
the molecular hydrogen column density derived from the 
1.3 mm continuum emission of Ward-Thompson et al. 
(1999), smoothed at a resolution of 22" (JV(H 2 )=1.1 x 10 23 
cm- 2 ): a;(H 2 D+) = 7V(H 2 D+)/A(H 2 ) - 4.4 x 10~ 10 . This 
value is only four times lower than the estimated electron 
abundance x(e) in the L1544 center (~ 2xl0" 9 , CWZ), 
suggesting that H 2 D + is one of the major ions in the gas 
phase (see Sect.^J. 

In case of subthermal emission, the line optical depth 
and the total column density increase. For example, if T ex 
= 6 K, t = 2.8 and the total H 2 D + column density be- 
comes ~ 1.2xl0 14 cm" 2 . This corresponds to ir(H 2 D + ) 
~ lxlO" 9 , so that x(H 2 D+) ~ x(e). We note that T cx = 
6 K is the minimum excitation temperature possible to 
have a solution of eq. JIJ in the case of an optically thick 
(r £ 0.5) line with T mb = 0.9 K. At the off-positions (20" 
or 2800 AU from the dust peak), where A(H 2 ) ~ 8xl0 22 
cm -2 , the derived H 2 D + abundance is 2.3xl0- 10 (if Tex 
= 7 K) and 4.0xl0" 10 (if T cx = 6 K), a factor of about 
two lower than toward the dust peak. 

We also carried out a radiative transfer 
analysis using the Monte Carlo program by 
Hogerheijde fc van der Tak (2000)f. The H 2 density 



profile was taken from lTafalla et al~l f2002k whereas dust 
temperatures are from lOalli et al.l f2002). Only thermal 
line broadening was included. Term energies, statistical 
weights and Einstein A coefficients were taken from 
the JPL catalog IjPickett et al. 1998jl . Besides collisional 
excitation, for which we used the scaled radiative rates 
proposed by Black et al. (1990) radiation from the 



cosmic microwave background and thermal radiation by 
local dust are taken into account, using grain opacities 
from Ossen kopf fc Henning| (1994). The abundance of 
H 2 D + was assumed to be constant within a radius ro 
and zero at r > rg. Model results are compared to line 
profiles at the center and offset positions to constrain 
the H 2 D + abundance and ro. The best- fit abundance is 
H 2 D + /H 2 = 10~ 9 , where calibration uncertainty allows a 
50% decrease. Models with r = 20-25 arcsec give good 
matches to the line fluxes at both the center and the 
offset positions. The optical depth at line center along the 
central pencil beam is 2.73. Here we limited our analysis 
to a static cloud. In a future paper (van der Tak et al., in 
prep.) we will explore the effects of including a velocity 
structure in the radiative transfer code. 
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http:/ /talisker. as. arizona.edu/~michiel/ratran. html 
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4. Chemistry 

The main result of this study is that H 2 D + has been de- 
tected for the first time in a pre-stellar core. The observed 
(±10 — In) line has a main beam temperature of ~ 1 
K, and the emission is concentrated in a region with ra- 
dius of ^2800 AU, roughly coincident with the so-called 
"flattened region" seen in the millimeter dust continuum 
emission, where the density is slowly changing with radius 
(~ r -14 ; see lAndre et al. 2f)0f)|l and/or the temperature 
is around 7 K IjEvans et al. 20011 lUalli et al.l 20021. 

The data have been analysed with a simple analytical 
model and with a Monte Carlo radiative transfer code. 
The derived H2D+ abundance (x(H2D + ) ~ 10~ 9 ) is con- 
sistent with the electron fraction estimated by CWZ in 
the center of L1544. This implies that H2D 4 " is the main 
molecular ion in the core center. The simple chemical code 
of CWZ, constrained by observations of several molecu- 
lar ions, predicts x(H2D + ) ~ 7xl0~ n in the center of 
L1544, about one order of magnitude less than that de- 
duced from the present data. A similar H2D" 1 " abundance 
(~ 6x 10~ n ) at the cloud center was also predicted by the 
detailed chemical model of lAikawa et al.l f2001'). coupled 
with the Larson-Penston dynamical evolution which best 
reproduced the L1544 observational results. 

CWZ found that to match the observed N 2 D+/N 2 H + 
(~0.24) and DCO+/HCO+ (-0.06) column density ratios 
toward the L1544 center one needs to (i) allow differential 
depletion of molecular species onto dust grains (in par- 
ticular, N2 has to be more volatile than CO, to maintain 
a large fraction of N2H + and N2D+ in the core center); 
and (ii) maintain a significant fraction of atomic oxygen 
in the gas phase, to limit the deuterium fractionation to 
the observed values (^ 0.2). As a by-product, Ii30 + was 
predicted to be abundant (x(H 3 + ) ~ 10~ 9 ) in the core 
center. 

The results of the present paper apparently contra- 
dict the CWZ conclusions, because the large H2D+ abun- 
dance can only be produced if all neutral species, includ- 
ing O and N 2 , are essentially frozen onto dust grains (see 
also |Bergin et al. 2002| ). In fact, only in this case can the 
H2D" 1 " abundance reach values close to the electron abun- 
dance. This can be seen with a simple chemical model, 
where the ingredients are H, H2, HD, H2D+, H^, CO, 
HCO + , DCO + and electrons. We neglect recombination 
onto dust grains, which is likely not to be an important 
process for the above molecular ions in the high den- 
sity core nucleus, where small grains and PAHs - signifi- 
cant carriers of negative charges ( |Lepp fc Dalgarno 1988| ) 
- are expected to be deposited onto bigger grains (e.g. 
|Ossenkopf fc Henn ing 1994). Following section 2 of Caselli 
(2002), assuming conditions appropriate for the central 
2500 AU of L1544 (kinetic temperature = 7 K, n(H 2 ) = 
10 6 cm" 3 , see lCalli et all 2002 s ). we have: 



where fcf is the forward rate coefficient of the proton- 
deuteron exchange reaction H^+HD H2D++H2 (see 
below) , and /d is defined as the ratio between the "canon- 
ical" CO abundance (9.5xl0~ 5 ; |Frerking et al. 1982| ) and 
the observed N (CO) /N(R 2 ). The H+ abundance can be 
estimated from the charge conservation equation, so that 
we obtain (with rate coefficients as in Caselli 2002): 



x(Kj 



(4) 



-Rdeut 



= x(H 2 D+) 
x(B+) 



2 10 1U fc f 



3107/d + x(e)/2 10" 



(3) 



{fux(e) + 4 10- 8 )(i?DEUT + 1) + 8 10- 8 

Equations 10) and directly furnish the H2D + and H3" 
abundances as a function of /d, once x(e) is known. Figure 
3 shows the result of this simple chemical analysis, with 
x(H 2 D + ) (thin curves) and -Rdeut (dashed curves) as a 
function of /d, assuming x(e) = 2xl0~ 9 (as roughly de- 
termined by CWZ in the L1544 center). We note that 
eqn. J3J) breaks down at large values of /d (<; 100), when 
dissociative recombination starts to dominate over proton 
transfer reactions. In the figure, two different values for fcf 
have been used: (i) the "standard" value fcf = 1.5xl0 -9 
cm 3 s" 1 (upper curves), and (ii) fcf = 3.5xl0 -10 cm 3 s -1 , 
recently measured bv lCerlich et al.l f 2002) (lower curves). 
The two data points in the figure are our present estimates 
of £(H2D + ), the arithmetic means of the abundances cal- 
culated with T ox =6 and 7 K (Sect. EJ): 7.2xl0" 10 and 
3.15xl0 -10 at the peak and off-peak positions, respec- 
tively. The corresponding /d values are from CWZ. The 
error bars indicate the uncertainties in the excitation tem- 
perature (—1 K; Sect. |2I) for estimating x(H2D + ), and cal- 
ibration errors affecting the determination of /d (about 
30%; CWZ). 

The observed H2D + abundances cannot be reproduced 
unless the CO depletion factor is —20 in the off position 
and ~40 at the peak (the "Monte Carlo abundance" indi- 
cates an even higher /d, —70), whereas the observed fu 
values are ~7 and ^9, at the same positions. Although our 
comprehension of the deuteration process in cold gas may 
be substantially incomplete, it is certain that measure- 
ments of /d are contaminated by the CO emission along 
the line of sight coming from regions with marginal CO 
depletion, so that the resultant /d's estimates are lower 
limits of the total amount of CO freeze out at the core 
center. Indeed, the chemical model of CWZ predicts an 
almost complete (<;99%) CO freeze out within the central 
2800 AU of L1544. Following the model results in Fig.|3 
the predicted /d values at the observed H2D+ abundances 
correspond to -Rdeut ~ 1 and 2 at the off and central po- 
sitions, respectively. 

If other neutrals reacting with H2D + and fpj", such 
as O and N2, are present in the gas phase, the predicted 
a;(H2D + ) abundance will drop as in the case of a low /d 
value (e.g. Caselli 2002). Thus, the present H2D 4 " observa- 
tions also require a heavy depletion of O (^99%) and N 2 
(~97%). A way to overcome the apparent inconsistency 
with CWZ (who need significant fractions of molecular 
nitrogen and atomic oxygen in the gas phase at r ~ 2800 
AU) and well reproduce both the moderate value (~ 0.2) 
of the 7V(N 2 D + )/7V(N 2 H + ) column density ratio as well 
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Fig. 3. (thin curves) x(H2D + ) vs. the depletion factor 
/d in the high density core nucleus of L1544 (see text). 
(dashed curves) Rbevt as a function of fu- Data points 
are the measured values at the peak and 20" off-peak po- 
sitions. Dotted lines indicate the -Rdeut and /d values 
predicted by the simple chemical model applied to the 
L1544 central 2800 AU (see Sectior0J), at the observed 
a;(H2D + ) values. Large H2D + abundances (> 10~ 9 ) are 
obtained only if /d > 70. 

as the large H2D+ abundance (which implies i?DEUT ~ 2) 
is to increase the depletion rate of O and N2 within the 
"flattened" region. This case requires central core densi- 
ties of about 10 7 cm~ 3 within a radius of 1000 AU, which 
is consistent with current mm continuum dust emission 
measurements if the central temperature drops to values 
of ^ 7 K (e.g. lEvans et al. 2001J1 . The conclusion is that 
present observations are compatible with H2D + being the 
most abundant molecular ion and with a total depletion 
of elements heavier than helium in the central 2800 
AU. This can also explain the lower H2D + abundance de- 
rived in t he direction of N GC1333-IRAS4 (x(R 2 B+) = 
3xl0~ 12 : IStark et al. 1999jl . where the young stellar ob- 
ject heats up the central zones, thus allowing the return 
of solid phase molecules back into the gas phase. A more 
comprehensive chemical model, coupled with a detailed 
radiative transfer code is however needed to better under- 
stand the whole data set available for L1544 and we are 
going to explore this in a future paper (van der Tak et al., 
in preparation). 
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